arXiv: 1503.08129vl [cond-mat.str-el] 27 Mar 2015 


Temperature and composition phase diagram in the iron-based ladder compounds 

Bai_^Cs^Fe2Se3 

Takafumi Hawai,^ Yusuke Nambu,^ Kenya Ohgushi,^’^ Fei Yasuyuki 

Hirata,^ Maxim Avdeev,^ Yostiiya Uwatoko,^ Yurina Sekine/ Hiroshi Fukazawa/ 

Jie Ma,^ Songxue Chi,^ Yutaka Ueda,^’^ Hideki Yoshizawa,^ and Taku J. Sato^ 

^Neutron Science Laboratory, Institute for Solid State Physics, 

University of Tokyo, Tokai, Iharaki 319-1106, Japan 

‘^Institute of Multidisciplinary Research for Advanced Materials, Tohoku University, Sendai 980-8577, Japan 
^Institute for Solid State Physics, University of Tokyo, Kashiwa, Chiba 277-8581, Japan 
^Department of Physics, Tohoku University, Sendai, Miyagi 980-8578, Japan 
^ Key Laboratory of Physics and Technology for Advanced Batteries (Ministry of Education), 

College of physics, Jilin University, Changchun, 130012, Peoples Republic of China 
^ Bragg Institute, Australian Nuclear Science and Technology Organization, 

Locked Bag 2001, Kirrawee DC NSW 2232, Australia 
^Quantum Beam Science Directorate, Japan Atomic Energy Agency, Tokai, Ibaraki 319-1195, Japan 
^Quantum Condensed Matter Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA 
^Toyota Physical and Chemical Research Institute, Nagakute, Aichi 480-1192, Japan 

(Dated: March 30 , 2015 ) 

We investigated the iron-based ladder compounds (Ba,Cs)Fe2Se3. Their parent compounds, 

BaFe2Se3 and CsFe2Se3, have different space groups, formal valences of Fe and magnetic structures. 

Electrical resistivity, specific heat, magnetic susceptibility. X-ray diffraction and powder neutron 
diffraction measurements were conducted to obtain temperature and composition phase diagram of 
this system. Block magnetism observed in BaFe2Se3 is drastically suppressed with Cs doping. In 
contrast, stripe magnetism observed in CsFe2Se3 is not so fragile against Ba doping. New type of 
magnetic structure appears in intermediate compositions, which is similar to stripe magnetism of 
CsFe2Se3, but inter-ladder spin configuration is different. Intermediate compounds show insulating 
behavior, nevertheless finite T-linear contribution in specific heat was obtained at low temperatures. 

BAGS numbers: 72.20.-i, 75.25.-j, 75.40.-s, 75.50.Ee 


I. INTRODUCTION 

Iron-based superconductors (SCs) have created consid¬ 
erable interests in condensed matter physics community 
since its discovery [l|. Despite intensive efforts, the su¬ 
perconducting mechanism is still controversial; spin fluc¬ 
tuation 0 , orbital fluctuation Q, or both of them are 
prime candidates for mechanism of the iron-based SCs. 
As these fluctuations are sensitive to crystal and mag¬ 
netic structure, it is desired to explore superconductivity 
in iron-based compounds with different crystal and mag¬ 
netic structures. 

To date iron-based SCs are known to have the two 
dimensional Fe square lattice with Fe atoms being tetra- 
hedrally coordinated by pnictogens or chalcogens. They 
show metallic behavior in electrical resistivity, and have 
antiferroniagnetic order 3, [HI]. The fill [ 1 , HI Bill 
and 122 [lO, [H systems have single-stripe magnetic 
structure and their magnetic moments are less than 1 
/iB- The 11 systems have double stripe magnetic struc¬ 
ture with ^ 2 /iB moment [mil- In contrast, the 245 
systems, A2Fe4Se5 {A = K, Rb, and Cs), are insulators 
and have two-dimensional ( 2 D) iron square lattice with 
periodic vacancy of Fe atoms [l 3 , [lH . They have block 
magnetism with magnetic moments ^ 3 /tb and have 
been considered to be the third type of magnetic struc¬ 
tures for parent compounds of iron-based SCs. 


To gain further insights into the interplay between 
crystal structure, magnetism and superconductivity, re¬ 
search on compounds with different dimensionality is im¬ 
portant. The spatial dimensionality should have effects 
on magnetic interaction, orbital selectivity and topology 
of Fermi-surface of Fe atoms. In this regard, it is notewor¬ 
thy that one-dimensional spin-ladder systems (La, Y, Sr, 
Ca)i4Cu2404i provide great insight into the understand- 
in^ of the two-dimensional copper-oxide superconductors 

Recently, iron-based ladder compounds AFe2Y3 {A = 
K, Rb, Cs or Ba, and X = chalcogen) dMi have been 
under intensive investigation | 2 Q|-[^. The AFe2Y3 com¬ 
pounds consist of FeX4 tetrahedra, which is similar to the 
2 D iron-based SCs, but Fe atoms form two-leg ladder lat¬ 
tice giving rise to the quasi-one dimensionality. They also 
show related magnetism with iron-based SCs, stripe and 
block magnetism. Unlike the most iron-based SCs, iron- 
based ladder compounds are insulators. KFe2Se3 and 
CsFe2Se3 [ 25 [ have stripe magnetism; magnetic moments 
are arranged to form ferromagnetic units along the rung 
direction, and the units stack antiferromagnetically along 
the leg direction. BaFe2Se3 BUI has block magnetism; 
magnetic moments are arranged to form Fe4 ferromag¬ 
netic units, and each units stack antiferromagnetically 
along the leg direction. Interestingly, in BaFe2Se3, Fe lat¬ 
tice distortion under Neel temperature (Tn) is reported 
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[^ . [ 2 ^ . This indicates that lattice and magnetism are 
strongly coupled. 

Recently magnetization measurement for the iron- 
based ladder compounds Bai_a^Ka^Fe 2 Se 3 have been re¬ 
ported [26|. BaFe 2 Se 3 has the Pnma space group, formal 
valence of Fe^+ and block magnetism with Tn ^ 255 K 
iSlii; on the other hand KFe 2 Se 3 has the Cmcm space 
group, formal valence of Fe^‘^+ and stripe magnetism 
with Tn ^ 200 K [26|. The previous study determined 
Tn of Bai_a^Ka,Fe 2 Se 3 from magnetization measurement 
and showed reduced Tn by K doping decreases [26| . How¬ 
ever, little is known on the variation of the crystal and 
magnetic structures in intermediate compositions. 

In this study, we focus on the mixed system of 
BaFe 2 Se 3 and CsFe 2 Se 3 . CsFe 2 Se 3 has the Cmcm space 
group, formal valence of Fe^‘^+, stripe magnetism similar 
to KFe 2 Se 3 with Tn ^ 175 K [25|. Since CsFe 2 Se 3 has 
lower Tn than KFe 2 Se 3 , we assumed that Tn was strongly 
suppressed by dilution. We examined Bai_a^Csa^Fe 2 Se 3 
through electrical resistivity, specific heat, magnetic sus¬ 
ceptibility, X-ray diffraction and powder neutron diffrac¬ 
tion measurements, and established composition and 
temperature phase diagram of Bai_a,Csa,Fe 2 Se 3 . The 
boundary of the two distinct space groups, Pnma and 
Cmcm, is between x = 0.05 and 0.15. Block mag¬ 
netism in BaFe 2 Se 3 is drastically suppressed by Cs dop¬ 
ing, whereas stripe magnetism in CsFe 2 Se 3 is relatively 
stable. In addition, new magnetic structure is observed 
for X = 0.5,0.55 and 0.65 at the lowest temperatures. 
Surprisingly at x = 0.25, no magnetic reflections were 
observed down to 7 K in the powder neutron diffraction. 


II. EXPERIMENTAL 

Polycrystalline samples of Bai_a,Csa,Fe 2 Se 3 {x = 
0.05,0.15,0.25,0.4,0.5,0.55,0.65 and 0.75) were synthe¬ 
sized from the mixture of BaFe 2 Se 3 and CsFe 2 Se 3 . The 
values of x were defined as nominal composition. The 
two parent compounds were synthesized from the stoi¬ 
chiometric ratio of the starting materials such as Ba : 
Fe : Se = 1 : 2 : 3 and Cs 2 Se : Fe : Se = 1 : 4 : 5 
[ 23 , [HI. They were mixed in a carbon crucible sealed 
into a quartz tube with 0.3 atm of Ar gas and reacted at 
900 °C for 24 hours. 

Single crystals of Bai_a,Csa,Fe 2 Se 3 {x = 
0.05,0.15,0.25,0.4,0.5,0.55,0.65 and 0.75) were grown 
from the melt. Stoichiometric mixtures of the two parent 
compounds were put in a carbon crucible and then sealed 
into a quartz tube with 0.3 atm of Ar gas. They were 
reacted at 1100 °C for 24 hours and cooled to 750 °C 
with a rate of 6 °C/h, and then cooled down to room 
temperature for 12 hours. Reflecting one dimensionality 
of the structure, the crystals preferentially grew along 
the leg direction. 

Quality of the samples was checked by the powder x- 
ray diffraction (XRD) with Cu-K^ radiation at room 
temperature. No detectable amount of impurities were 


TABLE I. Detailed diffractometer parameters. Samples are 
packed into vanadium cans. Temperature is controlled by 
using a Closed-cycle refrigerator (CCR). 


Diffractometer 

ECHIDNA 

WAND 

Wave length (A) 

2.4395(2) 

1.4827(1) 

26 range (deg.) 

2.75-163.875 

5.025-129.625 

20 step (deg.) 

0.125 

0.2 

Monochromator 

Ge(335) 

Ge(113) 

Refrigerator 

CCR (3-300 K) 

CCR (5-300 K) 


observed in the XRD measurements. Single-crystal XRD 
with Mo-Ka radiation was performed by using R-AXIS 
RAPID II (Rigaku) at room temperature to determine 
the space group of Bai_a,Csa,Fe 2 Se 3 . 

Electrical resistivity is measured using the four-probe 
method in the temperature (T) range between 30 and 
350 K. Specific heat measurements were performed using 
relaxation method between T = 1.8 and 300 K. These 
measurements were performed using the Physical Prop¬ 
erties Measurement System (Quantum Design, PPMS). 

DC susceptibility was measured between T = 2 and 300 
K under magnetic field up to = 5 T using a commercial 
SQUID magnetometer (Quantum Design, MPMS). 

Powder neutron diffraction was performed using the 
high-resolution ECHIDNA diffractometer at the Aus¬ 
tralian Nuclear Science and Technology Organization, 
and the Wide-Angle Neutron Diffractometer (WAND) 
at Oak Ridge National Laboratory. Experimental con¬ 
ditions are summarized in Table [H The Rietveld refine¬ 
ments were performed using the Fullprof suite [HI- The 
^stal structures were depicted by the VESTA software 


III. RESULTS 

A. CRYSTAL STRUCTURE 

Figure [T] shows neutron diffraction patterns for 
Bai_a,Csa,Fe 2 Se 3 together with Rietveld refinement re¬ 
sults at 300 K. 


Table HIl summarizes refined atomic positions, and other 
crystallographic data from the Rietveld refinement are 
shown in Table mil A little amount of carbon contami¬ 
nation was found for x = 0.65 and 0.75 samples, 0.55(1) 
% and 0.027(1) % weight ratio, respectively. In addition, 
an unindexed impurity phase was detected for x = 0.75 
sample, with the strongest peak which is roughly 20 times 
smaller than the main phase. Although their peak in¬ 
tensities are very weak, they could affect the refinement 
results. Hence, we excluded these 20 regions: 41.9-43.2° 
(carbon) for x = 0.65, and 24.5-25.3°, 47.8-48.7° (un¬ 
known impurity) and 42.2-42.9° (carbon) for x = 0.75 
during the refinements. 







3 



26 


FIG. 1. (color online). High-resolution powder neu¬ 
tron diffraction patterns of Bai_ccCsccFe 2 Se 3 at 300 K on 
ECHIDNA together with Rietveld analysis. The bottom blue 
lines give the difference between the observed (open circles) 
and calculated (line) intensities. 


The space group of Bai_a;Csa;Fe 2 Se 3 was deter¬ 
mined using single crystal XRD. We assumed that 
Bai_a^Csa^Fe 2 Se 3 belong to either Pnma or Cmcm^ and 
that the structure boundary should exist between x = 0 
and 1. It should be noted that the presence of 112Pnma 
nuclear reflection, which is forbidden in the Cmcm^ was 
used to determine the space group of Bai_a^Csa;Fe 2 Se 3 . 
We found that the 112Pnma reflection was present for 
X = 0.05, but absent for x > 0.15 in the rotation pho¬ 
tographs taken at the room temperature. It is reported 
that BaFe 2 Se 3 undergoes a phase transition from the 
Pnma to Cmcm structure at 660 K HH. Accordingly, it 
is natural to consider that the boundary of Pnma and 
Cmcm crosses 300 K between x = 0.05 and 0.15. 

Figure [Sfa) shows the crystal structure of 
Bai_a^Csa;Fe 2 Se 3 , and Fig. Hb) depicts one Fe 
ladder with both notation of Pnma and Cmcm. Figures 
me) and provide view being parallel to the leg 

direction of Cmcm and Pnma structure, respectively. 
The Fe ladders slightly incline in the Pnma structure. 
Since the definitions of a^b and c axes in Pnma and 
Cmcm are different, hereafter in this paper, the axes 
of Bai_a^Csa;Fe 2 Se 3 are defined as follows: the a axis is 


TABLE H. Rehned atomic positions at 300 K for 
Bai-ccCsa;Ee 2 Se 3 {x — 0.05, 0.15, 0.4, 0.5, 0.55, 0.65, and 
0.75). Note that x = 0.05 belongs to the Pnma space group. 
The others belong to Cmcm. 


Atom 

site 

X 

y 

Z 

t/iso (A^) 

X = 0.05 






Ba/Cs 

4c 

0.1798(8) 

1/4 

0.508(2) 

0.020(3) 

Fe 

8d 

0.4960(4) 

0.0008(10) 

0.3544(4) 

0.0120(9) 

Sel 

4c 

0.3632(6) 

1/4 

0.2265(8) 

0.022(3) 

Se2 

4c 

0.6281(6) 

1/4 

0.4959(11) 

0.0131(19) 

Se3 

4c 

0.3938(6) 

1/4 

0.8100(7) 

0.007(2) 

X = 0.15 






Ba/Cs 

4c 

1/2 

0.1788(8) 

1/4 

0.011(3) 

Fe 

8e 

0.3518(5) 

1/2 

0 

0.0223(11) 

Sel 

4c 

0 

0.1245(6) 

1/4 

0.011(2) 

Se2 

8^ 

0.2073(6) 

0.3804(5) 

1/4 

0.040(2) 

X = 0.4 






Ba/Cs 

4c 

1/2 

0.1725(8) 

1/4 

0.025(3) 

Fe 

8e 

0.3527(5) 

1/2 

0 

0.0185(10) 

Sel 

4c 

0 

0.1215(5) 

1/4 

0.0099(19) 

Se2 

8p 

0.2120(5) 

0.3818(5) 

1/4 

0.0429(19) 

X = 0.5 






Ba/Cs 

4c 

1/2 

0.1732(7) 

1/4 

0.023(3) 

Fe 

8e 

0.3540(5) 

1/2 

0 

0.0247(10) 

Sel 

4c 

0 

0.1195(5) 

1/4 

0.0163(19) 

Se2 

8fl 

0.2153(5) 

0.3833(5) 

1/4 

0.0467(19) 

X = 0.55 






Ba/Cs 

4c 

1/2 

0.1720(7) 

1/4 

0.028(3) 

Fe 

8e 

0.3539(4) 

1/2 

0 

0.0286(10) 

Sel 

4c 

0 

0.1209(5) 

1/4 

0.0149(17) 

Se2 

8^ 

0.2177(5) 

0.3837(5) 

1/4 

0.0522(19) 

X = 0.65 






Ba/Cs 

4c 

1/2 

0.1729(6) 

1/4 

0.028(2) 

Fe 

8e 

0.3552(4) 

1/2 

0 

0.0214(8) 

Sel 

4c 

0 

0.1200(4) 

1/4 

0.0175(15) 

Se2 

8fl 

0.2180(4) 

0.3840(4) 

1/4 

0.0404(14) 

X = 0.75 






Ba/Cs 

4c 

1/2 

0.1714(7) 

1/4 

0.031(2) 

Fe 

8e 

0.3563(3) 

1/2 

0 

0.0123(8) 

Sel 

4c 

0 

0.1205(4) 

1/4 

0.0103(16) 

Se2 

8p 

0.2208(4) 

0.3867(3) 

1/4 

0.0244(13) 


along the rung direction, the b axis is perpendicular to 
the two-leg ladder plane, and the c axis is along the leg 
direction, which is indeed the axis choice of the Cmcm 
space group. Figure [Sfa) shows x dependence of lattice 
constants a, b and c. Among the three lattice constants, 
the b value significantly deviates from the linearity in 
0 < X < 0.5. In particular for x = 0.15 roughly 2 
% deviation was found. It should be noted that this 
composition is close to the structural transition from the 
Pnma to the Cmcm. We also note that the decrease of 
the b value indicates reduction of inter-ladder distance, 
giving rise to enhanced three dimensionality. Figure [3l^b) 
shows that x dependence of local structure parameters at 
300 K. The definition of the local structure parameters 
are given in Fig. [2fb); the leg, rung, L indicate the Fe-Fe 
bond lengths; ai and 0^2 are the Se-Fe-Se bond angles. 
Local structures change continuously in x < 0.15 region. 
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TABLE III. Refined crystallographic data taken at 300 K. Chemical formula units Z is 4 for all x. Excluded regions contain 
peaks from carbon (20 = 41.9-43.2° for x = 0.65 and 42.2-42.9° for x = 0.75), and unknown impurity {20 = 24.5-25.3° and 

47.8-48.7 for X = 0.75). Rp = ll/obs 2/calc| / 2/obs5 Rwp — 'Wi j^obs 2/calc| / ^*2/obs] 


Formula 

X = 0.05 

(r = 0.15 

X = 0.4 

II 

0 

Oi 

(r = 0.55 

(r = 0.65 

X = 0.75 

Molar mass (g/mol) 

485.7 

485.2 

484.1 

483.7 

483.5 

483.0 

482.6 

Space group 

Pnma 

Cmcm 

Cmcm 

Cmcm 

Cmcm 

Cmcm 

Cmcm 

a (A) 

11.8648(3) 

9.2028(4) 

9.3242(3) 

9.3870(4) 

9.4083(3) 

9.4738(2) 

9.5885(2) 

b(A) 

5.47187(9) 

11.7954(5) 

11.8173(5) 

11.8152(5) 

11.8342(4) 

11.8412(3) 

11.8475(3) 

c {A) 

9.1775(2) 

5.5238(2) 

5.6091(2) 

5.6320(2) 

5.6413(1) 

5.6579(1) 

5.6785(1) 

Cell volume (A^) 

595.83(2) 

599.62^) 

618.045) 

624.645) 

628.095) 

634.7i5) 

645.075) 

F(OOO) 

191.65 

191.79 

192.14 

192.28 

192.35 

192.49 

192.63 

Calculated 

5.415 

5.375 

5.203 

5.144 

5.113 

5.055 

4.969 

density (g/cm^) 

Number of 

30 

23 

23 

23 

23 

23 

23 

parameters 

Excluded regions 

- 

- 

- 

- 

- 

41.9-43.2 

24.5-25.3; 

in 2 ^ (deg.) 

No. of observed 

199 

no 

112 

115 

115 

117 

42.2-42.9; 47.8-48.7 
118 

reflections 

Rp{%) 

14.2 

15.4 

15.0 

14.7 

15.2 

13.2 

12.7 

R^op{%) 

15.4 

17.4 

16.4 

15.5 

15.8 

14.0 

14.1 

Re.p{%) 

3.72 

3.78 

4.13 

4.42 

3.89 

3.96 

4.53 


17.1 

21.2 

15.8 

12.3 

16.5 

12.4 

9.75 


whereas weak anomaly could be seen around x = 0.05, 
which may be related to the structural phase transition. 
Figures [3l^c), [3fd) show temperature dependence of local 
structure parameters for x = 0.05. Although x = 0.05 
is Pnma^ being the same as BaFe 2 Se 3 , there is no 
enhancement of ladder distortion. Figures [3]^e), [3jf) 
show T dependence of local structure parameters for 
X = 0.75. Change in local structure was detected around 
T = 120 and 70 K, where two successive magnetic phase 
transitions were observed from the paramagnetism to 
stripe-II, and from the stripe-II to stripe-I, respectively 
(described in IIII Cp . 


B. BULK PROPERTIES 

Temperature dependence of electrical resistivity (p) 
along the leg direction is shown in Fig. SJ^a). No anomaly 
was seen in p in the whole measured temperature range. 
Intermediate compounds show much lower p than that of 
the parent compounds; at x = 0.55, p at 300 K shows the 
lowest value p = 1.89 x 10“^ U cm, which is 10^ times 
smaller than that of BaFe 2 Se 3 and 10 ^ smaller than that 
of CsFe 2 Se 3 . The temperature dependence of p shows 
one-dimensional variable-range-hopping (VRH) behavior 
in low-T region, given by p oc exp[(To/T)^/^] [Fig. lU^b)] 
[3ll-[^. It is known that {k^To)~^ of ID VRH is pro¬ 
portional to the localization length and the density of 
states at Fermi level. Figure [5] shows the (/cbTo)”^ ob¬ 
tained by fitting the data for Bai .^Cs.^Fe 2 Se 3 versus x. 
Although p of CsFe 2 Se 3 (x = 1 ) are deviated from 
the ID VRH behavior, we fit it with ID VRH to compare 
(/cbTo)”^ with those of Bai_a^Csa^Fe 2 Se 3 . Interestingly, 



FIG. 2. (color online). (a), (b) Crystal structures of 
Bai-ccCsccFe 2 Se 3 . A cuboid with solid lines indicates a crys¬ 
tallographic unit cell. Note that the definition of a, b and c 
axis is different from each other for the Pnma and Cmcm 
structures. Distinction for arrangement of ladder between (c) 
Cmcm and (d) Pnma space group. In Pnma, ladders are 
slightly tilted toward the a axis. 

the largest (/cbTq)”^ was observed for x = 0.65, which is 
closer to the strong insulator CsFe 2 Se 3 than BaFe 2 Se 3 . 

Figure |4l[c) shows specific heat {Cp) as a function of 
temperature. No anomaly was observed in all measured 
temperature and all x ranges. The CpjT versus 
plots at low temperature indicate T-linear contribution 
(7 term) extrapolated as T ^ 0 except for x = 0.15 
[Fig. lUd)]. This suggests that Bai_a,Csa,Fe 2 Se 3 have 
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T{K) _ T(K) 



FIG. 3. (color online), (a) Cs concentration (x) dependence 
of lattice constants at 300 K. The data for x = 0 and 0.05 are 
converted into the definition of Cmcm symmetry. The dashed 
lines represent linear fit to the lattice constants, (b) The x 
dependence of the Fe-Fe bond length along the leg and rung 
direction; Se-Fe-Se bond angle ai and 0 ^ 2 ; the ratio between 
the bond length along the leg and rung direction, and the 
diagonal Fe-Fe bond length (L) [see Fig. [2][b)]. (c), (d) T 
dependence of local structures for x = 0.05. Note that two 
distinct values for leg, leg/rung and ai are owing to Pnma 
space group (black square and rhomboid). The red circles 
show the average of the two values, (e), (f) T dependence of 
local structures for x = 0.75. Black lines indicate magnetic 
transition temperatures (see in text). 


finite 7 values. As discussed below, magnetization mea¬ 
surements indicate that the samples contain FeySeg. One 
can suspect that the obtained specific heat contain the 
contribution from the impurity. However, the amount 
of FcyScg is at most 1 % molar ratio (see below) and 
its 7 values are relatively small; 8 and 4.4 mJ/mol 
for 4c and 3c structures of FeySeg, respectively [sl, 1^ . 
Therefore the impurity effect is negligible. The 7 values 
obtained by linear fit for T < 3.3 K [solid lines in Fig. 
|4l[d)] are summarized in Fig. [5l We note that (/cpTo)”^ 
and 7 show a similar x dependence. 

Figure [ 6 ] shows the magnetic susceptibilities (x) de¬ 
fined as magnetization (M) over magnetic field (i7), that 
is, X = M/iT, for single crystals of Bai_a^Csa^Fe 2 Se 3 un¬ 
der H = 50000 Oe. The samples contain a tiny amount 



FIG. 4. (color online), (a) Temperature dependence of elec¬ 
trical resistivity of Bai-ccGsxFe 2 Se 3 along the leg direction, 
(b) One-dimensional variable range hopping type plot of (a). 
Solid lines indicate the results of fitting, (c) Temperature de¬ 
pendence of the specihc heat (Cp). Arrows show the magnetic 
transition temperature. Note that each data shifted by 100 
T/mol K for clarity, (d) Specific heat divided by tempera¬ 
ture against type plot. Solid lines indicate the results of 
fitting. Although Bai-ccGsxFe 2 Se 3 show an insulating behav¬ 
ior, 7 term (T-linear term) remains. 


of ferromagnetic impurity, and in this case, FeyScg is the 
most plausible one. Estimated molar ratio of FeySeg in 
Bai_a;Csa;Fe 2 Se 3 from hysteresis loop is less than 1 % for 
X = 0.05, and 0.5 % for other compositions [36|. We also 
measured the temperature dependence of magnetization 
for the annealed FeySeg between T = 2 and 300 K on 
field cool (EC) and zero-field cool (ZEC) measurements 
[inset of Eig. [ 6 ], and found that x decreases below 120 K 
on cooling. However x of Bai_a^Csa^Ee 2 Se 3 do not show 
such behavior, indicating that the negligible contribution 
from EeySeg impurity phase. The x of Bai_a^Csa;Ee 2 Se 3 
show paramagnetic increase as temperature decreases, 
and glassy behaviors were observed for 0.05 < x < 0.4 at 
low temperatures. The systematic decrease of the spin- 
glass transition temperature (Tj) implies that this be¬ 
havior is of intrinsic nature; if this behavior were due to 
ferromagnetic impurity phase, Tf would be constant. 


C. MAGNETIC STRUCTURES 

To elucidate variation of magnetic structures in 
Bai_a^Csa^Ee 2 Se 3 , we performed powder neutron diffrac¬ 
tion measurements. Eigures[7l [8]and[9]show the powder 
diffraction patterns of several temperatures from 3 to 300 
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FIG. 5. (color online). Cs concentration (a:) dependence of 
7 (left-axis) and (/cbTo)”^ (see text, right-axis). They show 
similar behavior and take the maximum value for x — 0.65. 
In background, the type of magnetism at low temperature is 
shown. 


K, in a selected scattering wave vector Q range. For 
X = 0.05,0.4,0.65 and 0.75 the diffraction patterns were 
collected using ECHIDNA, and for x = 0.15, 0.25, 0.5 and 
0.55 using WAND. 


1. 0<x< 0.15 

The diffraction patterns for x = 0.05 and 0.15 are 
shown in Fig. [71 For x = 0.05 magnetic reflection ap¬ 
pears at Q ^ 0.7 A“^, which corresponds to the main 
peak of block magnetism in BaFe 2 Se 3 appearing at the 
magnetic wave vector = (0.5, 0.5, 0.5) [^. This 
magnetic scattering appears below T* = 290(10) K, and 
grows with decreasing temperature. This peak is broad¬ 
ened beyond instrumental resolution even at the lowest 
temperature, suggesting that only 5 % doping eliminates 
the long ranged correlation of the block magnetism. Like¬ 
wise, for X = 0.15, only magnetic diffuse scattering was 
observed below T* = 215(15) K at similar position to 
the block magnetism. The correlation length is roughly 
estimated from the peak width as: 16(4) A for x = 0.05 
at 3 K, 9(3) A for X = 0.15 at 3.4 K. BaFesSeg {x = 0) 
have 19(5) A correlation length at 275 K [2^, thus these 
correlation lengths seem to decrease toward x = 0.25, 
where the magnetic orders are completely suppressed. 


0.4 < X < 1 

Figure [9] shows diffraction patterns for 0.4 < x < 1, 
showing stripe magnetism. We first focus on x = 0.75. 
On cooling, a magnetic diffuse scattering develops be¬ 
low 250 K at Q = 1.2 A“^. In contrast to the block 
magnetism, magnetic reflections become almost resolu¬ 
tion limited below 130(5) K. These magnetic reflections 



r(K) 

FIG. 6. (color online). Temperature (T) dependence of mag¬ 
netic susceptibility (x) for Bai-ccGsxFe 2 Se 3 along the a, 5 
and c axis in the Cmcm notation at the magnetic field H 
= 50000 Oe. The solid and open symbols indicate zero-field 
cooled (ZFG) and field cooled (FG) measurements, respec¬ 
tively. Note that each data are shifted for clarity. The ar¬ 
rows indicate the magnetic transition temperatures. The inset 
shows temperature dependence of y for the magnetic impurity 
FerSeg. 


change their positions below 80(5) K, which corresponds 
to additional magnetic transition [inset in Fig. fTOfa)]. 

All the magnetic peak positions below 80(5) K are well 
accounted for by the propagation vector Qm = (0.5, 0.5, 
0), which is the same as CsFe 2 Se 3 . To identify this mag¬ 
netic structure, we applied representation analysis and 
performed Rietveld analysis. Basis vectors (BVs) of the 
irreducible representations (irreps) of Qm were calculated 
using the SARA/i code The obtained BVs are listed 
in Table II in Ref. [25|. We sorted all BVs by comparing 
R-factors and found that T/^g, with spins aligned along 
the leg direction, has the best fit with Rp = 11.7 %; the 
second best is 12.8 % for '02- The Rietveld refinement 
results are shown in Fig. flOl a). The obtained magnetic 
structure is shown in Fig. fTTf b). We call this magnetic 
structure as stripe-I magnetism. This magnetic structure 
is completely the same as that of CsFe 2 Se 3 , implying the 
relatively stable stripe-I magnetism against the Ba sub¬ 
stitution. 

On the other hand, between 130(5) and 80(5) K mag¬ 
netic reflection positions cannot be accounted for by the 
propagation vector of BaFe 2 Se 3 , Qm = (0.5, 0.5, 0.5) nor 
of CsFe 2 Se 3 , Qm = (0.5, 0.5, 0). The peak positions are 
instead well accounted for by Qm = (0.5, 0, 0), however 
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FIG. 7. (color online). Powder neutron diffraction patterns 
for X = 0.05 and 0.15 taken at several temperature points. 
The data were collected at the listed temperatures using 
ECHIDNA for x = 0.05 and WAND for x = 0.15. 



FIG. 8. (color online). Powder neutron diffraction patterns 
for X = 0.25 taken at several temperature points. The data 
were collected at the listed temperatures using WAND. 


peak intensities of the magnetic reflections are too weak 
to perform Rietveld refinements. With decreasing x, we 
found that the new magnetic phase becomes the ground 
states in x = 0.5, 0.55 and 0.65 above 3 K and that the 
reflections are the most intense for x = 0.65. Therefore 
we analyzed the magnetic structure for x = 0.65. Figure 
me) shows a diffraction pattern for x = 0.65 taken at 3 
K and a result of Rietveld refinement. We performed rep¬ 
resentation analysis in the same way as x = 0.75 and ob¬ 
tained BVs summarized in Table HVl The sites of Fe atom 
in orbit 1 and orbit 2 cannot be exchanged by symmetry 
operations that leave the qm invariant. Therefore, we can 
choose different phase factors of moments for Fe atoms in 
orbit 1 and orbit 2 . Magnetic moment of the j-th Fe atom 
rrij {j = 1,2 G orbit 1 or j = 3,4 G orbit2) is given by 

the real part of ^^^orbitl ~ or 

"^i 6 orbit 2 = , where i 

means BV for j-th Fe atom belongs to orbit i {i = 
1 , 2 ), (pi is phase factor of orbit i, t is the lattice trans¬ 
lation vector and C is a real constant. We assumed that 
orbit 1 and orbit 2 are in the same irreps and that all 



FIG. 9. (color online). Powder neutron diffraction patterns 
for X = 0.4 taken at several temperature points, 0.5, 0.55, 
0.65, and 0.75. The data were collected at the listed temper¬ 
atures using EGHIDNA for x — 0.4, 0.65 and 0.75, and using 
WAND for X — 0.5 and 0.55. 


Fe atoms have the same moment size, because Fe atoms 
in paramagnetic state are in single site (see Table in|) . 
The best fit is a combination of with the phase factor 
(pi = 57 r /4 and (p 2 = 7 r/ 4 , where Rp = 12.0 % (the second 
best is 12.2 % for ^p^). 

The schematic magnetic structure is shown in Fig. 
me). We call this magnetic structure as stripe-II. Stripe- 
I and stripe-II have similar intra-ladder structure, where 
magnetic moments are arranged to form ferromagnetic 
units along the rung direction, and the units stack an- 
tiferromagnetically along the leg direction. In contrast, 
the spin directions [Fig. [TTl b) and fTH c)] and the inter¬ 
ladder relation [Fig. [TTT d) and [TTT e)] are different be¬ 
tween stripe-I and stripe-II. 

Knowing the magnetic structure of the x = 0.65 sam¬ 
ple, we tried to determine the magnetic structure for 
X = 0.75 at 90 K, where the same initial magnetic struc¬ 
ture is assumed. The best fits are obtained for ip^ and 7/^5 
with Rp = 11.7 %, however ipi have the close value with 
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TABLE IV. Basis vectors (BVs) of irreducible representations 
(irreps) for the space group Cmcm with the magnetic wave 
vector Qm — (0.5, 0, 0). Subscripts show the moment direc¬ 
tion. Columns for positions represent No.l: (a:,1/2, 0), No. 2: 
(a;,l/2,1/2), No. 3: (1-a;, 1/2,1/2) and No. 4: (1-a;, 1/2, 0). 


Irreps 

BV 

Orbit 1 

No. 1 No. 2 

Orbit 2 

No. 3 No. 4 

Ti 

'01 

2a 

-2a 

2a 

-2a 

Ts 

'02 

2a 

2a 

2a 

2a 

Vs 

03 

26 

26 

26 

26 


04 

2c 

-2c 

2c 

-2c 

r4 

05 

26 

-26 

26 

-26 


06 

2c 

2c 

2c 

2c 


Rp = 11.8 %. The X [Fig. [Hfb)] show little anisotropy 
for 0.5 < X < 0.75. Nevertheless it is hard to expect any 
anisotropic change without variations in significant struc¬ 
tural parameter for x = 0.65 and x = 0.75, we speculate 
that the spin orientation is parallel (t/^i) also at 90 K of 
the X = 0.75 sample. To confirm this orientation, single 
crystal neutron diffraction measurements is desired. 

For 0.4 < X < 0.75, magnetic diffuse scatterings were 
also observed at Q ^ 1.2 at higher temperatures. 
These scatterings evolve below ^ 200K [see rhomboids 
in Fig. ng. The positions are different from that of the 
X = 0.05 and 0.15 samples, where diffuse scatterings ap¬ 
pear at Q r\j 0.7 A- near the Bragg reflection position 
of the block magnetism. Therefore these diffuse scatter¬ 
ings would reflect stripe-type short-range correlation. No 
diffuse scattering was observed in CsFe 2 Se 3 . 


3. Magnetism of x — 0.25 

Interestingly, for x = 0.25, no magnetic reflection was 
observed down to 7 K as shown in Fig. [51 The y for 
X = 0.25 shows the spin-glass behavior, and it is con¬ 
sistent with this result. Note that other intermediate 
compounds (0.05 < x < 0.4) also show spin-glass behav¬ 
ior, however they have some magnetic scatterings (mostly 
diffuse) in powder neutron diffraction. 


4- Moment size of Bai-xCsxFe 2 Se 3 

Temperature dependence of estimated magnetic mo¬ 
ment sizes for several x are shown in Fig. [121 The mo¬ 
ment sizes of X = 0.65 and 0.75 are the results of Rietveld 
refinement. For x = 0.5 and 0.55 moment sizes were ob¬ 
tained from comparing between integrated intensity of 
the 021c'^cm nuclear peak and the ^21c'^cm magnetic 
peak. The x = 0.75 sample shows successive phase tran¬ 
sitions. Around 130 K the stripe-II magnetism appears 
and the magnetic moment grows to 0.31(3) /ip- In 85 - 
95 K, the stripe-I and stripe-II structures could coexist, 
because in this temperature range the moment size of 




FIG. 10. (color online). High-resolution powder neutron 
diffraction patterns for (a) x = 0.75 at 3 K, (b) x = 0.75 
at 90 K, and (c) x = 0.65 at 3 K using ECHIDNA with 
Rietveld analysis results (solid lines). The calculated posi¬ 
tions of nuclear and magnetic reflections are indicated (green 
ticks). The bottom line gives the difference. The insets show 
(a) change of positions of magnetic reflections from stripe-II 
to stripe-I magnetism and (c) temperature evolution of mag¬ 
netic reflections for x = 0.65. 


stripe-II becomes small. On further cooling, the stripe-I 
magnetism is finally stabilized, and the magnetic mo¬ 
ment reaches to 2.01(5) jxb at 3 K. The x = 0.5, 0.55 and 
0.65 samples, where stripe-II magnetism appears, show 
similar moment size to the x = 0.75 sample at the low¬ 
est temperature; 0.39(3) /ip for x = 0.5; 0.55(3) /ip for 
X = 0.55; and 0.60(2) /ip for x = 0.65. As x decreases to 
the region where magnetism is suppressed, the moment 
size decreases. 

In iron-based SCs, magnetic structure shows a close 
correlation with the ordered moment. Single-stripe mag¬ 
netism usually has an ordered moment smaller than 1 
/ip, double-stripe has ^ 2 /ip and block magnetism has 
a larger moment around 3 /ip. On the other hand for 
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FIG. 11. (color online). The schematic figures of magnetic 
structures for: (a) block magnetism (b) stripe-I magnetism, 
and (c) stripe-II magnetism. The inter-ladder relation of spins 
are shown in (d) stripe-I and (e) stripe-II magnetism. The 
+ and — sign indicate ferromagnetic and antiferromagnetic 
correlation in neighboring ladders, respectively. 


5. Magnetic phase diagram 

Figure [13] shows established temperature and compo¬ 
sition phase diagram of Bai_a^Csa^Fe 2 Se 3 . Interestingly, 
Tn is drastically suppressed by increasing x, and only 5 
% Cs doping (i.e. x = 0.05) destabilizes the long range 
order of the block magnetism. Almost simultaneously, 
space group changes from Pnma to the Cmcm between 
X = 0.05 and 0.15. Compared with the block magnetism, 
stripe-I and II magnetism can be seen in a large region 
of phase diagram and are relatively resilient the Cs dop¬ 
ing. For X = 0.75, successive magnetic phase transition 
from the paramagnetic state to stripe-II and from stripe- 
II to stripe-I was observed. With decreasing Cs concen¬ 
tration X, the stripe-II magnetism becomes ground state 
for 0.5 < X < 0.65 and Tn decreases gradually. The 
long range magnetic order is completely suppressed at 
X = 0.4. 

350f 
30d- 


250 r 

g200 

I- 150 
100 
50 




0.2 






oO' 


0.4 

X for Ba 




X ^ stripe-I magnetism - 

~ k Stripe-II 


0.6 

i-xCs/e2Se3 


0.8 


1.0 


iron-based ladder compounds, stripe-II magnetism has 
^ 0.5 /iB, stripe-I has ^ 2 /^b and block magnetism has 
< 3 /iB- The correlations are basically consistent be¬ 
tween the 2 D iron-based SCs and the ladder compounds. 
It may be noteworthy that the ladder compounds cover 
whole range of the moment size seen in iron-based SCs. 



FIG. 12. (color online). Temperature dependence of esti¬ 
mated magnetic moment. The data of x =0.65 and 0.75 were 
collected using EGHIDNA, and of x = 0.5 and 0.55 using 
WAND. 


FIG. 13. (color online). Magnetic phase diagram determined 
by the powder neutron diffraction. The data were collected 
down to 3 K (x = 0.05, 0.4, 0.65, and 0.75), 5K {x — 0.15, 0.5 
and 0.55), and 7 K (x = 0.25). Rhomboids indicate the tem¬ 
perature where diffuse scattering appears. Squares, upward 
triangles, and downward triangles indicate the magnetic phase 
transition temperature of the block, stripe-I, and stripe-II 
magnetism, respectively. 


IV. DISCUSSION 

We here discuss the interplay of crystal structure, bulk 
properties, and magnetism of Bai_a^Csa;Fe 2 Se 3 . 

The T-linear contribution (7 value) of Cp becomes 
large around x = 0.65 (Fig. |5j). Possible origin of fi¬ 
nite 7 values in Cp may be spin-glass state, magnons or 
electrons. 

Firstly, for the spin-glass possibility, we obtained the 
highest Tf ioT X = 0.15, and for 0.5 < x < 0.75 glassy 
behavior did not appear [Fig. | 6 ][b)]. This is not consis¬ 
tent with the behavior of 7 values, thus the possibility of 
spin-glass contribution can be excluded. 

Secondly, temperature dependence of Cp is known as 
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Cp (X for magnetic excitations associated with 

magnons, where d is spatial dimensionality, u = 1 for 
antiferromagnetic; = 2 for ferromagnetic ordered sys¬ 
tem. Thus T-linear contribution of the magnons could 
be due to ID antiferro-type or 2D ferro-type. The for¬ 
mer is plausible, and the latter should be unlikely due 
to antiferromagnetic structure of Bai_a,Csa,Fe 2 Se 3 . In 
stripe-II magnetism region, ordered magnetic moments 
are smaller than that of parent compounds and thus fluc¬ 
tuating moment may be significant. This part may con¬ 
tribute to the larger 7 value in the stripe-II region. How¬ 
ever, the parent compounds, which have larger moment 
size than intermediate compounds, do not show finite 7 
values. Therefore it would be difficult to conclude that 
the contribution of magnon is dominant. 

Finally, for the case of electronic specific heat, electrons 
may have density of states at Fermi level, but random po¬ 
tential prevents their itinerancy. In this possibility, the 
resistivity would follow the VRH model in low temper¬ 
atures, and the magnetic moment would reduce because 
of its itinerancy. This is consistent with the fact that 
similar fashion is seen in the 7 values and (/cbTo)”^ The 
small moments of stripe-II magnetism is also consistent 
with this scenario. These behaviors can be seen in in¬ 
sulating amorphous alloys, however theoretical model to 
account for 7 in insulating system is still controversial 
[ 3 ^. We suggest that this is the most plausible scenario 
among these possibilities. 

The obtained results indicate the close relation be¬ 
tween the crystal and magnetic structures. In the Pnma 
structure, the intra-ladder structure is distorted and 
rungs are slightly tilted from the parallel alignment along 
the rung direction [Fig. [2l[d)]. It is reported that 
BaFe 2 Se 3 shows enhancements of distortion in intra¬ 
ladder structure below Tn ^ 255 K, and this indicates 
magnetoelastic coupling On the other hand, at 

X = 0.05, where only magnetic diffuse scattering was ob¬ 
served, there is no drastic intra-ladder distortion in mea¬ 
sured temperature range [Fig. [3][c),[3l[d)]. In addition, at 
X = 0.15 where the crystal structure has the Cmcm space 
group, only weaker magnetic diffuse scattering than that 
oi X = 0.05 was observed. These facts suggest that no 
long range order of the block magnetism can exist in the 
Cmcm structure, and that the block magnetism may be 
sensitive to the distortion of intra- and inter-ladder atom 
arrangement. 

The stripe-I and II magnetism are more resilient than 
the block magnetism to the elemental substitution for 
0.4 < X < 0.75 region. No structure transition was ob¬ 
served in 0.4 < X < 0.75. The stability of the magnetic 
structure would be due to that of the crystal structure. 

For X = 0.75, the local structure changes slightly at the 
transition temperatures [Fig. [3][e), [3l[f)]. However the 
relation between the crystal and the magnetic structure 
is unclear at present. 

Next we focus on the x = 0.25 sample, where no dif¬ 
fuse or long range magnetic scattering observed by pow¬ 
der neutron diffraction measurements down to 7 K. In 


this case, moments would vanish or freeze without long 
range order nor short range correlation. From suscepti¬ 
bility measurements, the moments seem to remain finite 
[Fig. [ 6 ][b)]. Hence, there must be fluctuating moment (or 
glassy component) even at the lowest temperature. This 
compound has the competitions of block and stripe-II 
magnetism, as well as the Pnma and Cmcm structure. 
The suppression of magnetism for x = 0.25 would be due 
to these instabilities. 

We now compare obtained magnetic phase diagram 
(Fig. nsi) with theoretical calculations based on a five- 
orbital Hubbard model in a finite size Fe-Se two-leg- 
ladder [^. The calculations correspond to formal va¬ 
lences of Fe^+ (correspond to BaFe 2 Se 3 , or x = 0 ) and 
Fe2-25+ (Bao. 5 Cso. 5 Fe 2 Se 3 , 01 x = 0.5). They took ac¬ 
count of the lattice distortion seen in the Pnma space 
group as anisotropic hopping parameters, and they cal¬ 
culated the X = 0.5 sample as the Pnma (experimental 
results indicate x = 0.5 is Cmcm). Block and CX struc¬ 
ture [corre^onds to stripe-I and H magnetism, see Fig. 
1 in Ref. Q] are stable in a robust region of their phase 
diagrams for Fe2+ and Fe2-25+^ respectively. Except for 
the orientation of spins, their results are basically con¬ 
sistent with our findings. Interestingly, the stripe mag¬ 
netism can be stabilized in the Pnma space group. Hence 
the theory suggests that the stability of stripe-I and H 
magnetism are due to not only that of crystal structure 
but also that of electronic structure. Note that magnetic 
structures obtained by our study are realized in the re¬ 
alistic ratio J}i/U = 0.25 and U/W ^ 0.5 [see in Ref. 
[ 3 ^]. This implies Bai_a,Csa,Fe 2 Se 3 would be ’’interme¬ 
diate” coupling compounds 

Pressure-induced superconductivity was reported in 
cuprate ladder compounds Sri 4 _a,Caa,Cu 2404 i [iol-li^. 
This compound has CU 2 O 3 combined two-leg-ladder 
structure, and is insulator. Their resistivity becomes 
lower with increasing pressure and show superconductiv¬ 
ity at 3 GPa. Moreover, the suppression of magnetism 
provides superconductivity in 2 D iron-based SCs, and 
we discovered that the magnetism of Bao. 75 Cso. 25 Fe 2 Se 3 
{x = 0.25) is completely suppressed. Therefore applying 
pressure for this composition may be an intriguing way to 
pursue superconductivity in future. Pressure-effect study 
is now in progress. 


V. CONCLUSION 

We performed a comprehensive study on the interplay 
of crystal structure, bulk properties, and magnetism for 
Bai_a,Csa,Fe 2 Se 3 by electrical resistivity, specific heat, 
magnetic susceptibility. X-ray diffraction and powder 
neutron diffraction measurements. The space group 
changes from Pnma to Cmcm between x = 0.05 and 
0.15, where block magnetism is drastically suppressed 
concomitantly. New type of magnetic structure, stripe-II, 
appeared where T-linear term of Cp is finite. For x = 0.75 
we observed successive phase transitions from stripe-II to 
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stripe-I with a slight structure distortion. Typical mag¬ 
netic moment sizes in Bai_a;Csa;Fe 2 Se 3 was determined 
as: stripe-II magnetism ^ 0.5/iB; stripe-I ^ 2/iB; and 
block magnetism < S/zb- They cover the moment size 
range seen in the all iron-based SCs: single-stripe mag¬ 
netism < I/tb; double-stripe ^2/iB; and block magnetism 
^ 3/iB- Notably, there is no diffuse nor resolution-limited 
magnetic scattering in powder neutron diffraction pat¬ 
tern down to 7 K for X = 0.25, where competition of the 
block and the stripe-II magnetism and of the Pnma and 
Cmcm structure is expected. 


The work at IMRAM is also supported by Cooperative 
Research Program of ’’Network Joint Research Center 
for Materials and Devices”, and ”Nano-Macro Materials, 
Devices and System Research Alliance”. The research 
at Oak Ridge National Laboratorys High Flux Isotope 
Reactor was sponsored by the Scientific User Facilities, 
Office of Basic Energy Sciences, U.S. Department of En¬ 
ergy. 


ACKNOWLEDGMENTS 

The work at ISSP and IMRAM is supported by KAK- 
ENHI (26800175, 23340097, 23244068, and 24224009). 


[1] Y. Kamihara, T. Watanabe, M. Hirano, and H. Hosono, 
J. Am. Chem. Soc. 130, 3296 (2008) 

[2] K. Kuroki, S. Onari, R. Arita, H. Usui, 

Y. Tanaka, H. Kontani, and H. Aoki, 

Phys. Rev. Lett. 101, 087004 (2008), 

[3] H. Kontani and S. Onari, 

Phys. Rev. Lett. 104, 157001 (2010) 

[4] M. D. Lumsden and A. D. Christianson, 
J. Phys.: Condens. Matter 22, 203203 (2010) 

[5] P. Dai, J. Hu, and E. Dagotto, Nat. Phys. 8, 709 (2012) 

[6] H. Luetkens, H.-H. Klauss, M. Kraken, F. Litterst, 
T. Dellmann, R. Klingeler, C. Hess, R. Khasanov, A. Am¬ 
ato, C. Baines, et al., Nat. Mater. 8, 305 (2009). 

[7] J. Zhao, Q. Huang, C. de La Cruz, S. Li, J. Lynn, 
Y. Chen, M. Green, G. Ghen, G. Li, Z. Li, et aL, Nat. 
Mater. 7, 953 (2008). 

[8] G. Ghu, F. Ghen, M. Gooch, A. Guloy, B. Lorenz, 
B. Lv, K. Sasmal, Z. Tang, J. Tapp, and Y. Xue, 
Physica G 469, 326 (2009), superconductivity in Iron- 
Pnictides. 

[9] S. Li, G. de la Gruz, Q. Huang, G. F. Ghen, T.- 

L. Xia, J. L. Luo, N. L. Wang, and P. Dai, 
Phys. Rev. B 80, 020504 (2009) 

[10] Q. Huang, Y. Qiu, W. Bao, M. A. Green, J. W. Lynn, 
Y. G. Gasparovic, T. Wu, G. Wu, and X. H. Ghen, 
Phys. Rev. Lett. 101, 257003 (2008) 

[11] J. Zhao, W. Ratcliff, J. W. Lynn, G. F. Ghen, 
J. L. Luo, N. L. Wang, J. Hu, and P. Dai, 
Phys. Rev. B 78, 140504 (2008) 

[12] W. Bao, Y. Qiu, Q. Huang, M. A. Green, P. Zajdel, 

M. R. Fitzsimmons, M. Zhernenkov, S. Ghang, M. Fang, 
B. Qian, E. K. Vehstedt, J. Yang, H. M. Pham, L. Spinu, 
and Z. Q. Mao, Phys. Rev. Lett. 102, 247001 (2009) 

[13] O. J. Lipscombe, G. F. Ghen, G. Fang, T. G. 
Perring, D. L. Abernathy, A. D. Ghristianson, 
T. Egami, N. Wang, J. Hu, and P. Dai, 
Phys. Rev. Lett. 106, 057004 (2011) 

[14] J. Guo, S. Jin, G. Wang, S. Wang, K. Zhu, T. Zhou, 
M. He, and X. Ghen, Phys. Rev. B 82, 180520 (2010) 

[15] F. Ye, S. Ghi, W. Bao, X. F. Wang, J. J. Ying, 
X. H. Ghen, H. D. Wang, G. H. Dong, and M. Fang, 


Phys. Rev. Lett. 107, 137003 (2011) 

[16] T. Vuletic, B. Korin-Hamzic, T. Ivek, S. Tomic, 

B. Gorshunov, M. Dressel, and J. Akimitsu, 

Phys. Rep. 428, 169 (2006) 

[17] H. Hong and H. Steinfink, 

J. Solid State Ghem. 5, 93 (1972) 

[18] K. Klepp, W. Sparlinek, and H. Boiler, 

J. Alloys Gompd. 238, 1 (1996) 

[19] R. H. Mitchell, K. G. Ross, and E. G. Potter, 
J. Solid State Ghem. 177, 1867 (2004) 

[20] Y. Nambu, K. Ohgushi, S. Suzuki, F. Du, 

M. Avdeev, Y. Uwatoko, K. Munakata, 

H. Fukazawa, S. Ghi, Y. Ueda, and T. J. Sato, 

Phys. Rev. B 85, 064413 (2012) 

[21] V. Svitlyk, D. Ghernyshov, E. Pomjakushina, 

A. Krzton-Maziopa, K. Gonder, V. Pom- 
jakushin, R. Pottgen, and V. Dmitriev, 

J. Phys.: Gondens. Matter 25, 315403 (2013) 

[22] B. Saparov, S. Galder, B. Sipos, H. Gao, S. Ghi, D. J. 
Singh, A. D. Ghristianson, M. D. Lumsden, and A. S. 
Sefat, Phys. Rev. B 84, 245132 (2011) 

[23] J. M. Garon, J. R. Neilson, D. G. Miller, A. Llobet, and 
T. M. McQueen, Phys. Rev. B 84, 180409 (2011) 

[24] H. Lei, H. Ryu, A. 1. Frenkel, and G. Petrovic, 
Phys. Rev. B 84, 214511 (2011) 

[25] F. Du, K. Ohgushi, Y. Nambu, T. Kawakami, M. Avdeev, 
Y. Hirata, Y. Watanabe, T. J. Sato, and Y. Ueda, 
Phys. Rev. B 85, 214436 (2012) 

[26] J. M. Garon, J. R. Neilson, D. G. Miller, 

K. Arpino, A. Llobet, and T. M. McQueen, 
Phys. Rev. B 85, 180405 (2012) 

[27] J.-K. Bao, G.-M. Feng, Y.-K. Luo, H. Jiang, Y.-L. Sun, 
W.-H. Jiao, G.-Y. Shen, Z.-A. Xu, and G.-H. Gao, 
J. Phys.: Gondens. Matter 26, 026002 (2014) 

[28] E. Dagotto, Rev. Mod. Phys. 85, 849 (2013) 

[29] J. Rodriguez-Garvajal, Physica B 192, 55 (1993). 

[30] K. Momma and F. Izumi, 

J. Appl. Grystallogr. 44, 1272 (2011) 

[31] N. F. Mott, Conduction in Non-Crystalline Materials, 
2nd ed. (Glarendon Press Oxford, 1993). 

[32] E. M. Hamilton, Philos. Mag. 26, 1043 (1972) 



12 


[33] B. L Shklovskii and A. L. Efros, Electronic Properties of 
Doped Semiconductors (Springer-Verlag Berlin, 1984). 

[34] A. Okazaki, J. Phys. Soc. Jpn. 16, 1162 (1961) 

[35] K. Shimada, T. Mizokawa, K. Mamiya, T. Saitoh, A. Fu¬ 
jimori, K. Ono, A. Kakizaki, T. Ishii, M. Shirai, and 
T. Kamimura, Phys. Rev. B 57, 8845 (1998) 

[36] To estimate the effect of ferrimagnetism of FerSeg, we 
took hysteresis loops of the samples and pelleted FerSeg 
between H = -10000 and 10000 Oe at room tempera¬ 
ture, and compared with that of Bai_xCsa;Fe 2 Se 3 . We 
had annealed FerSeg using the same sequences as for 
Bai_ccCsa;Fe 2 Se 3 before measurements and measure it, 
since FerSeg has two superstructures 4c and 3c [53] . 

[37] A. Wills, Physica B 276-278, 680 (2000) 


[38] A. Y. Rogatchev and U. Mizutani, 
Phys. Rev. B 61, 15550 (2000) 

[39] Q. Luo, A. Nicholson, J. Rincon, S. Liang, J. Riera, 

G. Alvarez, L. Wang, W. Ku, G. D. Samolyuk, A. Moreo, 
and E. Dagotto, Phys. Rev. B 87, 024404 (2013) 

[40] M. Uehara, T. Nagata, J. Akimitsu, 

H. Takahashi, N. Mori, and K. Kinoshita, 
J. Phys. Soc. Jpn. 65, 2764 (1996), 

[41] T. Nagata, M. Uehara, J. Goto, J. Akimitsu, N. Mo- 
toyama, H. Eisaki, S. Uchida, H. Takahashi, T. Nakan- 
ishi, and N. Mori, Phys. Rev. Lett. 81, 1090 (1998), 

[42] M. Isobe, T. Ohta, M. Onoda, F. Izumi, S. Nakano, J. Q. 
Li, Y. Matsui, E. Takayama-Muromachi, T. Matsumoto, 
and H. Hayakawa, Phys. Rev. B 57, 613 (1998) 


